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We have observed 3C 279 in a 7-ray flaring state in November 2008. We construct quasi- 
simultaneous spectral energy distributions (SEDs) of the source for the flaring period of 2008 and 
during a quiescent period in May 2010. Data have been compiled from observations with Fermi, 
Swift, RXTE, the VLBA, and various ground-based optical and radio telescopes. The objective 
is to comprehend the correspondence between the flux and polarization variations observed during 
these two time periods by carrying out a detailed spectral analyses of 3C 279 in the internal shock 
scenario, and gain insights into the role of intrinsic parameters and interplay of synchrotron and 
inverse Compton radiation processes responsible for the two states. As a first step, we have used a 
multi-slice time-dependent leptonic jet model, in the framework of the internal shock scenario, with 
radiation feedback to simulate the SED of 3C 279 observed in an optical high state in early 2006. 
We have used physical jet parameters obtained from the VLBA monitoring to guide our modeling 
efforts. We briefly discuss the effects of intrinsic parameters and various radiation processes in 
producing the resultant SED. 



INTRODUCTION 



Blazars are well known for their variability and 
power of polarized radiation across a wide range of 
the electromagnetic spectrum [ll| [24], HH. In some 
cases, the flux can vary on timescales as short as 
an hour or less [see e.g., 0, [22j]. Blazars exhibit a 
doubly-peaked spectral energy distribution (SED), in 
which the low-energy component could extend from 
radio through UV or X-rays while the high-energy 
component extends from X-rays to 7-rays. The low- 
frequency component of the SED is almost certainly 
due to synchrotron emission from nonthermal, ultra- 
relativistic electrons. The high-frequency component, 
on the other hand, is a result of inverse Compton scat- 
tering of seed photons by the same ultra-relativistic 
electrons producing synchrotron emission (in a lep- 
tonic jet scenario). In this case, the seed photons 
could be the synchrotron photons produced within the 
jet (synchrotron self Compton, SSC) [2l|,[3lj], and/or 
external photons entering the jet from outside (the EC 
process) [e.g.,[H, and references therein]. The spectral 
variability patterns and SEDs are key ingredients in 
determining the acceleration of particles and the time- 
dependent interplay of various radiation mechanisms 
responsible for the observed emission. 

Another defining characteristic of blazars is the high 
degree of linear polarization at optical wavelengths. 
Many bright 7-ray blazars that arc in the Fermi- L AT 
Bright 7-Ray Source List [2J have shown spectral and 
linear polarization variability flEl . [LSI . l30j . Linear po- 
larization at mm, IR, and optical wavelengths tends to 
exhibit similar position angles and sometimes correla- 
tion across these wavebands, often with some time de- 
lay [HI, [l^, H3, [H, 28] . Such correspondence between 



the variation in polarization and flux across a wide 
range of the electromagnetic spectrum, combined with 
VLBI imaging, can be used to identify the location of 
variable emission at all wavebands and shed light on 
the physical processes responsible for the variability 

The blazar 3C 279, located at a redshift of 0.538 
[TTI |. is one of the most prominent and well-studied 
blazars. This is due to its highly variable nature 
(change in magnitude Am ~ 5 at optical bands) at 
all wavelengths and high optical polarization up to 
45.5% observed in the U band |32|. Intensive multi- 
wavelength campaigns [see, e.g.. [la Hil |27| and theo- 
retical efforts [e.g.,0, Qi| have led to some important 
conclusions about the physical properties of 3C 279. 
Chatterjee et al. (2008) showed that the flux vari- 
ability in 3C 279 has been found to be significantly 
correlated at X-rays, optical R band, and 14.5 GHz 
wave bands, which also suggests that nearly all X- 
rays are produced in the jet. The X-ray flux has also 
been associated with superluminal knots, as suggested 
by correlation with the flux of the core region in the 
43 GHz VLBA images [l2T |. Nevertheless, the nature 
and origin of its high-energy emission and the rela- 
tionship of its behavior to the physical aspects of the 
jet remain elusive [I3,|27j- In addition, the correspon- 
dence between high- and low-energy emission is also 
not very well understood. 

Here, we aim to understand the physical state of 
3C 279 at different flux levels and look for a correspon- 
dence between the flux and polarization variation ob- 
served during the flaring state of November 2008 and 
quiescent state of May 2010. This can be achieved 
by carrying out a detailed spectral analyses of 3C 279 
under the internal shock scenario, and understanding 
the role of intrinsic parameters and interplay of syn- 
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chrotron and inverse Compton radiation processes in 
shaping the corresponding spectra of the two states. 

As a first step toward comparing the physical state 
of 3C 279 at different flux levels, we use the multi-slice 
time-dependent leptonic jet model of Joshi & Bottcher 
(2011) (in the framework of internal shock scenario) 
with radiation feedback to simulate the SED of 3C 279 
corresponding to the optical high state of early 2006. 
The broadband emission of 3C 279 for this state indi- 
cates suppressed external Compton emission, which 
makes it an ideal candidate for simulation using a 
synchrotron-SSC model. We use physical jet parame- 
ters obtained from the VLBA monitoring to guide our 
modeling efforts and discuss the role of various intrin- 
sic parameters and radiation processes in producing 
the resultant SED. 

We briefly describe the model of Joshi & Bottcher 
(2011) in iJTTJ We discuss our findings about the con- 
nection between the flux and polarization variation 
observed in 3C 279 during flaring and quiescent pe- 
riods of November 2008 and May 2010, respectively, 
in £11111 We discuss our first results from this study in 
i ]III A[ We summarize and give a brief description of 
future work in ^IVI 
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FIG. 1: Schematic of the emission region with RS trav- 
elling into the inner shell and FS moving into the outer 
shell. The bulk Lorentz factors (BLFs) of the two shells 
are Fj & F , respectively, such that Fi > r o . The primed 
quantities refer to the comoving frame and unprimed re- 
fer to the lab (AGN) frame. The comoving pressures p' {s 
and Prs of the shocked fluids across the CD are considered 
equal. A rs and Af s are the widths of the inner and outer 
shells, respectively, after the collision. The quantities have 
been obtained by using hydro - dynamic and relativistic 
jump conditions across the shocks and CD @, Hj. 



II. INTERNAL SHOCK MODEL 

The mode of acceleration of particles to highly rel- 
ativistic energies and its location in the jet is not yet 
completely understood. A colliding shell or an inter- 
nal shock model offers a way to gain insight into the 
physics of particle acceleration. Under this model, 
the central engine (black hole + accretion disk) is as- 
sumed to spew out shells of plasma with different ve- 
locity, mass, and energy. An inner shell, ejected at a 
later time, and travelling with a higher speed catches 
up to a slower moving outer shell that is ejected at 
an earlier time. The subsequent collision results in 
an emission region, as shown in Figure [T] with two 
internal shocks (reverse (RS) and forward (FS)) sep- 
arated by a contact discontinuity (CD) and traveling 
in opposite directions to each other in the frame of 
the shocked fluid. As the shocks propagate through 
the emission region, they convert the ordered bulk ki- 
netic energy of the plasma into the magnetic field en- 
ergy and random kinetic energy of the particles. The 
highly accelerated particles then radiate and produce 
the emission observed from the jet. The treatment of 
shell collision and shock propagation is hydrodynamic 
and relativistic in nature [1, Hj| • 

The model follows the evolution of the electron 
and photon populations inside the emission region in 
a time-dependent manner. The equations governing 
this evolution are, 



dn e (^,t) 
dt 



d_ 



+ Qe(7,t) - 



(1) 



and 



dn ph (e,t) 
dt 



"ftph^em i^-i ^ph^abs 



n p h{e,t) 



ph,esc 



(2) 

Here, (d-j/dt)i oss is the particle (referred to as elec- 
trons in the rest of the text) energy loss rate due to 
synchrotron and inverse Compton losses. The quan- 
tity Qe{"f,t) is the sum of external injection of elec- 
trons and intrinsic 7 — 7 pair production rate, and 
£ e ,esc is the electron escape time scale. The quantities 
n p h, C m(e, t) and n p h,abs(e;^) are the subsequent pho- 
ton emission and absorption rates, and £ p h,csc is the 
photon escape timescale (see ^11 A[) . The radiative en- 
ergy loss rates and photon emissivities are calculated 
using the time-dependent radiation transfer code of 
Joshi & Bottcher (2011). 

The model reproduces the NIR to 7-ray emission 
from blazars as it follows the evolution of the emis- 
sion region out to sub-pc scales only. Thus, it simu- 
lates the early phase of 7-ray production. During this 
time, the radiative cooling is strongly dominant over 
adiabatic cooling and the emission region is highly 
optically thick at GHz radio frequencies. Hence, the 
simulated radio flux is well below that of the actual 
radio data. 
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FIG. 2: Schematic of the radiative transfer in between 
the zones using the appropriate photon escape probability 
function. The unprimed n ph ^ fwd / back values represent the 
photon densities in the forward and backward direction, 
respectively, in the comoving frame of the emission region. 
The rest of the unprimed quantities refer to the lab frame. 

A. Multi-zone Radiation Transfer Method 

We calculate the resultant spectrum from the col- 
lision, for a cylindrical emission region, in a time- 
dependent manner. We consider the inhomogeneity 
in the photon and electron energy density through- 
out the emission region by dividing the region into 
multiple zones and providing a fraction of the photon 
flux from each zone to its adjacent zones in the two 
directions, as shown in Figured) 

The fraction of photon flux from a zone or, in other 
words, the photon escape rate of a zone in a particular 
direction is calculated using the photon density of that 
zone and the probability of escape, P, for a photon 
from that zone in that direction (forward, backward, 
or sideways), according to the equation 



dn p h,fwd/back/sidc(e, ^) _ n p h(e, fi) 



dt 



P 



ph,csc 



fwd/back / side 



where t p \ l csc is the volume-and-angle-averaged photon 
escape timescale for a cylindrical region that is cal- 
culated semi-analytically in the model. We use this 
framework to calculate the radiation transfer within 
each zone and in between the zones [2631 ] . 



III. THE BLAZAR 3C 279 

Analytical studies of the SEDs of 3C 279, col- 
lected during the lifetime of CGRO, indicate that the 
SSC mechanism is responsible for producing X-rays 
[Tol l29ll . Intensive multiwavelength campaigns [see 
e -g- IS 13 H3 an d theoretical efforts [see e.g. 0, [13 
have provided insights on the physical properties of 
3C 279. But many aspects of the radiation mechanism 
responsible for high-energy emission still remain un- 
clear. Also, no consistent trends of cross-correlations 
have been found between optical, X-ray, and 7-rays 
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FIG. 3: Observed lightcurves of 3C 279 from 08/09/2008 - 
04/16/2011 in 7- and X-ray energy bands. The top panel 
corresponds to the data compiled from RXTE while the 
bottom panel refers to the Fermi data. Long-term vari- 
ability can be seen in the X-rays & 7-rays. 

[23| although the source, in early 2009, exhibited co- 
incidence of a 7-ray flare with that in the optical and 
a dramatic change of optical polarization angle during 
that time [l[ . Such behavior poses a serious challenge 
to the single-zone leptonic jet model approach where 
the electron population is assumed to be homogeneous 
and is responsible for producing synchrotron emission 
as well as high-energy emission through inverse Comp- 
ton scattering [l3j . 

In order to improve our insight into the behavior 
of 3C 279, we aim to understand the evolution of 
the source at different flux levels by simulating the 
corresponding SEDs, using a multi-zone leptonic jet 
model, and looking for connections between the evo- 
lution of spectral states and the behavior of the cor- 
responding polarization at optical and radio wave- 
lengths. Figure [3] shows the longterm 7- and X-ray 
lightcurves (LCs) of 3C 279 from August 2008 till 
April 2011. From these LCs, we have extracted two 
1-month time periods: Fl (11/08/2008- 12/08/2008) 
and Ql (05/22/2010 - 06/26/2010), corresponding to 
a 7-ray flaring and a quiescent period, as shown in 
Figures H&E1 

As can be seen from Fig. [3l long-term variability 
exists both in the X-rays & 7-rays. On the other hand, 
in case of Fl (see Fig. [4] more complicated features 
exists for such shorter time-periods that need to be 
considered in detail. While no correlation is apparent 
between the flux and polarization behavior at optical 
and radio wavelengths for period Fl, the VLBA data 
on 06/14/2010 (period Ql) resulted in a degree of po- 
larization p — 4.7 ± 0.8% and a radio polarization 
angle EVP A = 78.3 ± 4.9°, which agrees very well 
with the corresponding values of p (8.6 ± 0.3%) and 
EVPA (75.8 ± 0.9°) at optical R band. This suggests 
a common source of origin for the two. 

The SEDs of 3C 279 corresponding to states Fl 
and Ql are shown in Fig. HI The data have been 
compiled from observations with Fermi, Swift, VLBA, 
and various ground-based optical and radio telescopes. 
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FIG. 5: Multi-frequency lightcurves for the period Ql. ri- 
band fluxes correlate with the optical polarization. 

The optical and near-IR data were dereddened us- 
ing the Galactic extinction coefficients of Schlcgcl ct 
al. (1998) and the magnitudes were converted to 
fluxes using the zero-point normalizations of Bessel 
et al. (1998). The shape of the optical and near-IR 
spectra indicates that the synchrotron component for 
both states peaks in the near-IR regime with the vF v 
value of Fl higher than that of Ql by almost an or- 
der of magnitude. The X-ray flux of Fl remains in 
a low state coinciding with that of state Ql. The 
shape of the optical spectrum suggests that the syn- 
chrotron component extends into the UV regime with 
a turnover to the high-energy component taking place 
at ~ 10 16 Hz. The 7-ray spectrum for both states in- 
dicates a higher total energy output in the high-energy 
component compared to the low-energy one, which is 
typical for the class of flat spectrum radio quasars Q . 



A. First Results 

The blazar 3C 279 was observed in an exten- 
sive multi- frequency campaign in January 2006 (l3j . 
The object went into an optical flaring state on 
01/15/2006. The interesting feature of this state is 
that the X-ray flux remained in a low state during 
that period, while coinciding with the historical X- 
ray low-state of June 2003 [l3j]. Almost a month 
later, in Februrary 2006, it was detected by Major At- 
mospheric Gamma-ray Imaging Cherenkov Telescope 
(MAGIC) in a very high energy (VHE) 7-ray state [3j. 

As a first step toward simulating the SEDs of 3C 279 




v[Hz] 

FIG. 6: The observed SEDs of 3C 279 corresponding to 
periods Fl and Ql. 

in order to understand the evolution of spectral states 
and compare the flaring states of the source in Novem- 
ber 2008 and January 2006, we have used the time- 
dependent leptonic jet model of Joshi & Bottcher 
(2011) that calculates the synchrotron and SSC radi- 
ation to reproduce the observed SED of 3C 279 corre- 
sponding to the optical high state of 01/15/2006. Fig- 
ures [7|& [5] show the instantaneous and time-integrated 
simulated SEDs of 3C 279, respectively, for this day. 
The model independent paramete rs [9j estimated us- 
ing the SED, VLB A observations [2J|, and variability 
on 1-day timescale were used to develop an initial set 
of input parameters: 



S 
R 
D 

Tmin 
7max 

p 

^obs 



15.5 
2.5 x 



10 



16 



0.82 4 A G 
9.0 x 10 2 
2.2 x 10 4 
4.3 

2.1° ±1.1° 



(4) 



Here, 6 is the Doppler boosting factor. The symbol 
R is the radius of the emission region, B is the comov- 
ing magnetic field of the region, and £b is the equipar- 
tition parameter between the magnetic field and the 
electron energy density, assumed to be equal to 1 here. 
The quantities 7 m i n & 7 max refer to the low and high 
energy cutoffs of the electron energy distribution. The 
spectral index of the electron population, at the time 
of injection, is given by q and p is the equilibrium 
spectral index obtained using the optical synchrotron 
spectrum, which is given by p = q + 1 for strongly 
cooled electrons. The quantity, # b s is the observing 
angle inferred from VLBA observations. 

The entire emission region has been divided into 100 
zones (50 in the forward and 50 in the reverse emission 
region) to analyze the observed SED of the source. 
The initial set of parameters was modified to repro- 
duce the state of 3C 279 as observed on 01/15/2006. 
Various simulations were carried out by varying in- 
put parameters one at a time to obtain a reasonable 
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FIG. 7: Simulated instantaneous SEDs of 3C 279 for 
01/15/2006 showing the time-dependent evolution of the 
optical high state. Each solid curve is representative of the 
spectrum corresponding to a particular instant of time. 

fit to the observed SED of 3C 279. A fit is con- 
sidered successful if it passes through the observed 
data points, corresponding to 01/15/2006, without 
over-producing the X-ray flux. An observer's line of 
sight angle, 9 b s = 1-6°, and a cross-sectional radius, 
R = 6.85 x 10 16 cm, for the cylindrical emission region, 
have been chosen to obtain the current best fit. Ta- 
ble U lists the rest of the parameters used for obtaining 
the current fit of the observed SED. These parameters 
result in r s h ~ 18.6 for the entire emission region, 
B « 1.3 G and 7 ma x ~ 1-2 x 10 5 for both the forward 
and reverse emission regions, and 7 m i n . f s ~ 1-4 x 10 3 
& 7mm, rs ~ 1.7 x 10 3 for the forward and reverse 
emission regions, respectively. 

As can be seen from Figures [7] & [51 the observed 
SED for 01/15/2006 shows a high-energy bump that 
is indicative of a dominant SSC component and a sup- 
pressed EC component. The time-integrated simu- 
lated SED of the successful model passes very close 
to the IR and optical data points, indicating that the 
synchrotron component is responsible for the lower 
energy bump of the SED. The spectral upturn takes 
place in the soft X-rays at > 0.1 keV due to the pres- 
ence of the SSC component in the simulation. The 
lower-energy part of the SSC component reproduces 
the X-ray data quite well, suggesting the dominance 
of the SSC component in producing this part of the 
high-energy bump. As can be seen in Figure[51 there is 
a break in the low-energy part of the X-ray spectrum 
of January 2006 at ~ 10 19 Hz, which could be due 
to the presence of a low-energy cutoff of the electron 
distribution at ultrarelativistic energies fi~3l ]. 



IV. DISCUSSION AND FUTURE WORK 

We have presented first results of our study on the 
behavior of 3C 279 at different spectral states, and the 
connection between the flaring and quiescent states 
and the behavior of its polarization at optical and ra- 




FIG. 8: Simulated time-integrated SED, averaged over 1 
day, of 3C 279 for 01/15/2006. 

TABLE I: Model parameters used to reproduce the state 
of 3C 279 as observed on 01/15/2006. 
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L w : luminosity of the injected electron population in the 

blob, ri, : BLFs of the inner and outer shells before 
collision, q: particle spectral index, e e : ratio of electron 
and shock energy density, eb'- ratio of magnetic field and 
shock energy density, £ e : fraction of accelerated electrons 



dio wavelengths. 

In order to carry this out, we have extracted a 7- 
ray flaring (Fl) and a quiescent (Ql) state from the 
long-term LCs of 3C 279 shown in Fig. [3] In case of 
Fl, a clear cross-correlation cannot be seen at optical 
frequencies over a shorter time-period such as a month 
(see Fig. [3} and some complicated features exist that 
need to be considered in detail. In case of Ql, there is 
an anti-correlation between 7- and X-ray, with X-rays 
mildly flaring during the 7-ray quiescent state, but the 
R-band flux shows a strong relationship to the optical 
polarization (see Fig. [5]). On the other hand, while 
no correlation can be seen between the flux and po- 
larization behavior at optical and radio wavelengths 
for the period Fl, the VLBA data during the Ql pe- 
riod agrees very well with the corresponding values 
of p and EVP A at optical R band. This suggests a 
common source of origin for the two during this state. 

In order to understand the behavior of the source 
at different spectral states, it is necessary to simulate 
the corresponding SEDs of the source. To that end, 
we have simulated the optical high state of 3C 279 
observed in January 2006 using the time-dependent 
multi-zone leptonic jet model of [26| . The model, in 
its current form, calculates radiation resulting from 
synchrotron and SSC processes. The spectral state 
of 3C 279 observed in January 2006 makes a good 
candidate for this model, as the 7-ray emission is sup- 
pressed in this case and the X-ray emission is generally 
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associated with SSC scattering pjj HH [29| . 

Our current model of the SED results in a high value 
of L w compared with that of Collmar et al. (2010). 
This is reasonable, as the internal shock model is ex- 
pected to have low acceleration efficiency, and thus 
would need a higher kinetic luminosity to achieve the 
desired acceleration. The current fit results in a value 
of magnetic field that agrees reasonably well with the 
model independent parameter estimate and that of 
Collmar et al. (2010). The value of electron energy 
index, q, used to obtain the current fit is almost equal 
to the model independent parameter p and is com- 
paratively higher than that of Collmar et al. (2010). 
This indicates that the synchrtron spectrum in the 
current fit does not come from a cooled population of 
electrons. The value of q needs further adjustment in 
order to obtain a satisfactory fit. With a successful fit, 
we hope to gain insights on the acceleration efficiency, 
the relative speeds with which the two shocks are mov- 
ing away from each other (in the frame of the shocked 
fluid), and the mass of each plasma shell needed to 
reproduce an optical high state of 3C 279 such as that 
of 01/15/2006. 

As part of our future work, the external Compton 
component due to photons entering the jet from the 
BLR and dusty torus is being incorporated in the ex- 



isting model of [26] to reproduce more fully the SEDs 
of blazars, especially for flat spectrum radio quasars 
(Joshi et al. 2012, in prep.), and to analyze a wider 
variety of blazar spectral states. We plan to use the 
extended version of the model to simulate the SEDs of 
3C 279 for a typical period of 10 days corresponding 
to Fl and Ql and compare the results with that of 
the January 2006 optical high state of 3C 279. 

Further, we plan to incorporate the effects of mag- 
netic field orientation, as inferred from polarization 
monitoring programs, on the resultant spectral vari- 
ability and SEDs of blazars. This would aid us in the 
study of intrinsic parameter differences between var- 
ious blazar subclasses that arise from the orientation 
of the magnetic field in the jet. 
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